TITLES 

TYPES 

DBSCRIPTIONs 



CAPACITYs 



METHOD OF USE: 



UHiVERSiry OF illhiois 

DIGITAL COMPUTER 

LIBRARY ROUTIH E K l6 - 263 

Modified Multiple Regression Analysis (DOI Only) 
Entire program 

This routine, which is a revision of K ik, will 
calciilate and print any or all of the following? 

(1) Zero-order correlations, 

(2) Covariances, 

(5) Means and standard deviations, 

(k) Multiple correlations with associated regression 
weights, both standardized and imstandardized, 

(5) The standard errors of the regression weights and 
the standard error of estimate. 

This analysis can be done on any subset of the variables 
on the data tape, and in addition, the user has a 
choice of five transformations that can be made on 
the variables. 

This routine differs from K ik in that the raw scores 
are stored on the drum. For repeated use, this routine 
is preferred as it is necessary to read the data once 
only. 

s 2 san^jle size 

ps total number of variables on tape 
n: subset of variables used in the analysis 
s(p + 1) < 8157; p < 50; n < 22; s > n + 1 

Stops 

1. Master tape 2^012 

2. Parameter tape. Type N 2^113 
3o Data tape 2iK)12 
There are two types of parameter tapes: 

(1) Type N directs the computer to read data from tape, 
store information on the drum, and compute. 

(2) Type J directs the computer to read from the drum 
and compute. 

Any number of Type J tapes can follow step 5 above 
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and each vlil Stop on 21*012. The data on the dnni 
can he replaced at any 2^1012 stop hy reading a Type 
H tape followed by a data tape as in steps 2 and 5 
above. 
PARAMETER TAPES, TYPE Type H and Type J parameter tapes are alike in eveiy 
K 4XD T3CPE J: respect ^cept for the final terminating symbol which 

for Type N is an N and for Type J is a J. 
The first section of each parameter tape must contain 
six parameters, s, n, t,i , k, and J. These are 
punched on tape with terminating symbols as follows: 

sSnHfF^LkKjJ 
(The terminating symbols, S and K, can also be 
punched as - and +). 

(1) "s" is the sample size or the number of observations. 

(2) "h" is the nuBttber of non-eliminated variables; i.e., 
the subset selected by the parameter tape for analysis • 

(3) "f" is the number of decimal places for the zero-order 
correlations. To suppress this print-out, set fa 0. 

(I4.) «£" is the number of decimal places for the covariances 
To suppress this print-out, set jf =0. 

(5) "k" is the number of places for the multiple 
correlations, the standardized regression weights and 
their standard errors. To suppress this print-out, 
set k = 0. 

(6) "j" is the number of places for the standard errors 
of estimate, the unstandardized regression weights^and 
their standard errors. To suppress these, set J = 0. 
The second section of each parameter tape consists of a 
se<iuence of I's and O's pimched on tape in the 
following way: 

00101011 space 
These digits correspondi to the n non-eliminated variables. 
The I's specify that these variables are to be 
considered successively as dependent, with the other 
(n - 1) variables In the subset as independent. 



- 5 - 



variables . la the exai^le alaore, the third, the 
fifth, the seventh, and the eighth vould he tr^ted 
in turn as dependent variables. The sequence Is 
terainated by any fifth-hole character such as a 
space, a carriage return, or a delay. It ia not 
necessary to punch O's after the final "1" in the 
sequence. As an example, if only the first variable 
is to be considered as dependent, then punch only 
a "1" followed by a space. 

The third section of each parameter tape is the 
identification. This can be a name, a date, a 
problem description, or any coanbination of teletype 
characters whatever except single fifth-hole delays. 
This portion of the tape is reproduced exactly as a 
heading for the results. It is terminated by a 
single hole delay. An exai^ae of such a heading 
follovs: 

Vl/59 

PROBLEM BO. 17 Single fifth-hole delay 
If no heading is desired, punch tvo fifth-hole delays. 
The fourth section of each parameter tape specifies 
the transfonaations to be made on the variables and 
whether a variable is to be eliminated or not. This 
section consists of a set of p signed integers where 
these integers correspond to the numbers on the data 
tapeo If X is the variable, the meaning of the 
integers is shown in the chart below s 
Integer Meaning 

+0 Eliminate x 

+1 X is unchanged 

+2 X is transf oimed to x 

+5 X is transformed to x 

•HI- X is transformed to the square root of x 

+5 X is transformed to '^'^J^ 

(Notes X must be positive for the +k and +5 
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DATA TAPE' 



SCALING OF THE DATA* 



transformations, and x must be greater than 

10~ for the +5 transformation . ) 
The set of signed integers must be terminated by an 
M for a Type N parameter tape and by an NJ for a 
Type J tape. An exanrple may be as follows; 

+1+0+0+1+1+1+2+2+0+1NN 
The p variables are punched on tape as signed fractions 
with an N terminating symbol after each of the s 
observations. If an F is punched instead of an N, 
the computer will stop on 5012K. An additional 
section of the data tape can be inserted in the 
reader! ^y raising the black switch, the problem 
is continued. 

The data are read from tape as signed fractions. 
Locating the decimal point in the unstandardized 
numbers can be confusing imless all of the variables 
are scaled by the same constant. 
If the scaling constant for all x^ is lO' , then 
the scaling in the results is as follows* 
(l) The means, standard deviations, and the standard 



error of estimate are scaled by 10 



-z 



.-2z 



(2) The covariances are scaled by 10 """, and 
• '(3) Tlie remainder of the results are unaffected^ 
MATHEMATICAL FOEMULAS: Let the jth observation of the ith variable be x, . 

j. 1, 2, ^. , s 



k; 1, 2, 



n 



s 
The iieaz^s x. = E 


a 




The Covariance: 


9 

E 


8 


The CarrelBtlan Coefflelent: 


=ik 






^ n — c — 
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The Standardized Regression Weight: 



"h- 



11 



ik 



vhere r Is the (l>k) element of the inverse of the 
correlation matrix . 

The Standard Error of the Standardized Regression Weights: 



Ik 



11 kk , ikv2 
(s - n) (T^^f 



The Multiple Correlation Coefficient; 



% = 



1 - 



11 



The Standard Error of Estimate: 



^1 = 



=11 (1 - <) 



The Unstandardized Regression Weights: 



Ik 



'11 



'kk 



^!k 



The Standard Error of the Unstandardized Regression 
Weights: 



% 



Ik \ 



11 



'kk 



^&. 



THE RESULTS: 



The form of print-out of the results can "be 
illustrated by an example. For a problem vlth 
four variables, the ^rlnt-out could appear as 
f ollovs . The mean and standard deviations are 
always printed to nine places. 
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?2sr^ 



CknnrelfttlQits. 




HflBIS • 



*0,16493l^^ 





"+-21^150000 

•»-*5'^96O0O00 
■I- .105000000 



•«• .00159^ 
+•0559559»^ 

+#0065(KllT 
•i>«008f{)l09 

♦ .0080l}^5 
-.OO55W0O 
♦•00089155 
+•00332500 

♦•117765966 

♦ ^36550079 
+ .051970115 
+ •057662813 



m^ ■' — Covarlffttetf 



^ Standardi drriAtioni 



Variable 



CoBstaat tezm. 



DURATION: 



+ #Ol^38l3563 
+.928217019 



-Multiple correlatiQii 



■^ -l^OOOOOOOO +.000000000 

-0.35097^0 +.07629333^ 

42 •92808*^3^ +.332283964 

-0.72196344 1-^.2778155451 

— »^ 1+0.160047881 "l '^'^^' 

tftistajidai^lzed Standard 
fiegresslon Weights Errors 



-1.00000000 

-0.66480972 

+1.29216327 

[-0^35pi4ll 



+.000000000 
+.153246236 
+.146636873 



\f'l^' 



standardized 
s slots Weights 



Standard 
Errors 



The chart below may be helpful in making time estimates. 
These are average times for typical problems , If the 
number of variables, p_, on the data tape, exceeds the 
number selected for the analysis, n, by a large amount, 
the reading time should be increased proportionately. 
AloO the number of digits on input and output should 
be taken into account. 

TIME IN sec::m)£ 



Variables 


Read Time (5 Digits) 


Calculation 


Printing Time ("^ Digits) 


n 


Type N 


Type J 


Invert R 


R 


Gov. 


Reg. Weights 


5 


.25 s 


.10 s 


k 


3 


r 


4 


10 


.'2 s 


.22 s 


10 


10 


10 


8 


15 


.88 s 


M s 


18 


20 


20 


12 


20 


1.50 s 


.70 s 


35 


35 


35 


17 
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HOTE If 



NOTE 2: 



NOTE 5S 



If the master tape stops on FF from location 578^ 

the sum check has failed indicating the tape has 

been read incorrectly. 

If the nuniber of variables, p, in any rov of the 

data tape does not agree with the number of integers 

on the parameter tape, the computer will stop on 

FF from location 12J. 

If a drum error has been made, the coaiputer will 

stop on FF from location 06n. 



RT: 1/15/60 



DAT E April 6, 19^9 



ROUTINE K Ik BY G.H. Golub 



REVISION BY V. West 



WRITE-UP BY K.W. Dickman 



APPROVED B Y %^ A^AJ<j^/C>Ui^ 



Igr 
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COMPUaJKR CffiDERS AHD HOPESs ' This routine is a revision of K ik. Changes 

in K 14 are shovn helow. 



LOCATICar ORDER 



NOTES 



Abs. Rel 



I 



Subroutine B Ccarputes means and cross-products divided hy s 



205 5 50L 267L Remove stop 

215 8 2637IF OOF Transfer to read final term, symbol 



• • • o • 



Routine H 



Read in data and perform transformations 



275 5O5OSJ 50L Same 

276 1 26sk 26297F Transfer to new tape stop 

290 15 JHF 22219F Bypass old tape stop 

291 16 261IL OOI6L Same 



297 


22 


4013F 


L012S5 


298 


23 


30298F 


26366F 


299 


Zh 


L521Si^ 


lo6f 


500 


25 


40F 


L3F 


501 


26 


36302F 


FFI8F 


302 


27 


26372F 


oof 


305 


28 


LJ12F 


5637L 






OO366K 




366 





L5328F 


0020F 


%1 


1 


k6k-L 


4615L 


368 


2 


JO5OSJ 


502L 


569 


3 


26S9 


OOIO8OS7 


370 


k 


OOF 


26299F 


371 


5 


4125L 


2618L 


372 


6 


L5528F 


IA5S5 


373 


7 


Fif3L 


403L 


37^ 


8 


F525L 


iK)25L 


375 


9 


lo4f 


36275F 



Stop on F^ transfer to store on dnmi 

Stop on FF if p variables on data tape do not 
agree with integers on parameter tape 

Same 



Store data on drum with Y 1 



Transfer to test number of variables read 



Increase drum address by p + 1 
Count rows 
Tests S rows? 



/ 
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LOCAXIQI 


QRXSES 




lOTES 




Ata. 


Bel. 


.««.•••••« 




^.™— ^ — ..^.-— ———*.— Mi«^————. 




376 


10 


PPIOL 


L512L 






577 


11 


ii0276r 


???76f 






378 


12 


2613L 


OOF 






579 


13 


5050SJ 


5015L 


Read frcM dnni using T 1 




380 


Ik 


2689 


001080S7 






581 


15 


OOF 


F532aF 






582 


16 


IA5S5 


lJu^l 


Xncrenae drum address 




385 


17 


HOI^L 


26277P 


Tranaf or« row 




381* 


18 


I523L 


kC21*:>¥ 






585 


19 


QlkT 


L02^L 


Read tendjiatlxig Bymbol 




586 


20 


362IL 


2437 


Stops WTd of paraaeter tape 




587 


21 


L532aF 


602OF 






388 


?? 


ij6l5L 


2210L 


Set nuiiber of varlahles to he read 




389 


23 


26^ 


OOIF 






390 


2k 


OOF 


0013F 








Sum Check Change 








L50F6 


89289 










2in8!r 









